Space-time codes with controllable ML decoding 
complexity for any number of transmit antennas 



Naresh Sharma Pavan R. Pinnamraju Constantinos B. Papadias 

Tata Institute of Fundamental Research Indian Institute of Technology Athens Information Technology 

Mumbai, India New Delhi, India Athens, Greece 

nareshs@ieee.org pavanramesh.iitd@gmail.com papadias@ait.edu.gr 



Abstract — We construct a class of linear space-time block codes 
for any number of transmit antennas that have controllable 
ML decoding complexity with a maximum rate of 1 symbol 
per channel use. The decoding complexity for Al transmit 
antennas can be varied from ML decoding of 2^ log2 M l _1 symbols 
together to single symbol ML decoding. For ML decoding of 
2 rio g2 Mi-n (n = 1,2, ••■) symbols together, a diversity of 
min(M, 2^'° S2 M l-™+ 1 ) can b e achieved. Numerical results show 
that the performance of the constructed code when 2 ^ log2 1 
symbols are decoded together is quite close to the performance 
of ideal rate-1 orthogonal codes (that are non-existent for more 
than 2 transmit antennas). 

I. Introduction 

Multiple antenna systems have been of great interest in 
recent times, because of their ability to support higher data 
rates at the same bandwidth and noise conditions; see e.g. 
[1],[2], [3], [13] and references therein. While orthogonal 
designs offer full diversity with single symbol ML decoding, 
they don't have rate 1 for more than 2 transmit antennas. 

The loss of rate has been addressed by the use of quasi- 
orthogonal codes that make the groups of symbols orthogonal 
where each group has more than one symbol in general [7], 
[8], [10], [12]. A fully orthogonal code would have just one 
symbol per group. Because of this relaxation of constraints, 
these codes achieve higher code rates that were hitherto not 
possible with orthogonal codes. It was shown in [9], [11], [14], 
[15] that performance of above quasi-orthogonal codes can be 
improved with constellation rotation. 

Codes for any number of transmit antennas were presented 
in [12]. In this paper, we construct that a new class of space- 
time codes with a maximum code rate of 1, that are inspired 
from the codes in [12], that have a useful property that the 
ML decoding is controllable. On one extreme, one can design 
rate 1 codes that have single symbol ML decoding offering 
diversity of 2, and on the other, one can have codes offering 
full diversity with ML decoding of M/2 symbols together. 

It is, however, shown for the constructed codes that for rate 
one codes with single symbol ML decoding, full-diversity is 
impossible and for codes that require more than one symbols 
to be decoded together for ML symbols decoding, it is indeed 
possible to have full-diversity. 

We use the following notation throughout the paper: *, T 
and f denote the conjugate, transpose and conjugate transpose 
respectively of a matrix or a vector; I M and M are M x 



M identity and null matrices respectively; || A \\p, det(A) 
and Tr(A) denote Frobenius norm, determinant and Trace of 
matrix A respectively; C denotes the complex number field; 
CAf (0,1) denotes a circularly symmetric complex Gaussian 
variable with zero mean and unit variance. 

II. System Model and Design Criterion 

Consider a system of M transmit and N receive antennas. 
For the ease of presentation, in this paper, we will assume 
that M is a power of 2. The case of M not being a power of 
2 can be treated easily as in [12] by constructing a code of 
size 2 r iog 2 M l and deleting columns suitably chosen to have 
the code matrix of size 2^2 M l x M. 

The statistically independent modulated information sym- 
bols are taken P at a time denoted by c = (ci, • • • , cp) T . This 
information vector is pre-coded (i.e. multiplied) by a M x P 
matrix denoted by 1Z. Let s = (si, • • ■ , s M ) T and 

s = Kc (1) 

with E{|s,;| 2 } = 1, i = 1, • • • , M. As we shall soon see, the 
choice of 1Z is central to the construction of codes, s is the 
input to a linear space-time block code that outputs a M x M 
matrix Gp[s], where 

M 

Gm[s] = (C m s m + D m s* m ) , (2) 

m— 1 

where C m , D m , m = 1, • • • , M, are MxM complex matrices, 
which completely specify the code. This code is transmitted 
in M channel uses and the average code rate is hence P/M 
symbols per channel use. For a quasi-static fading channel, 
the received signal is given by 

X[s} = ^G M [s]H + V, (3) 

where X and V are the M x N received and noise matrices, 
and H is the M X N complex channel matrix that is assumed 
to be constant over M channel uses and varies independently 
over the next M channel uses and so on. The entries of H 
and V are assumed to be mutually independent and CJ\f(0, 1), 
and p is the average SNR per received antenna. We assume 
that channel is perfectly known at the receiver but is unknown 
at the transmitter. 



It has been shown in [1] by examining the pair- wise 
probability of error between two distinct information vectors 
(say u, v e C p ), that for full-diversity, in quasi-static fading 
channels, G M [JZ(u — v)]Gm[K(u — v)] should have a rank 
of M. For square code matrices, the above criterion could be 
modified to yield 



min det{G M [^(u - v)]} ^ 

U.V.U^V 



(4) 



III. Iterative construction of space-time codes 

The main difference between these codes and those in [12] 
is the choice of 1Z that will allow us to vary the ML decoding 
complexity and construct full-diversity codes with decoding 
of a pair of symbols. 

Let us define two disjoint partition vectors that are function 
of the vector s (whose length will be clear from the context) 
denoted by _4m,i(s) and Am, 2 (s)- These partition vectors 
have same length as s and have the same symbols as s in 
indices they possess and zeros in other indices. If we denote 
the first and last M elements of a 2M x 1 vector s by 
Sm,i and sm,2 respectively, then these partitions are iteratively 
constructed as 

-4.2M,i(s) = -4m,i(sm,i) + Am,2(sm,2) (5) 
-42M,2(s) = -4m,2(sm,i) + -4m,i(sm,2), (6) 

and the code is iteratively constructed for the zth partition as 



G 2M [A 2M ,i(s)] 



Gm [Am,i{sm,1 )} Gm [■A-M'i ( S M,2 )] 

-G m [-4m,;( s m,2)] G M [A M ,i( s *MA)] 



(7) 

where i = 2, if i = 1 and is 1 otherwise, and hence by using 
linearity, we have 



G 2 m\s] = 



Gm[sm,i] Gm[sm,2] 
-Gm[s* M2 \ Gm[s*mi\ 



(8) 



where Gi[s] = si V s e C 1 , Ai : \ = s\, and A 2 ,i is a null 
set. 

A. Receiver Processing 

We give a practical decoding algorithm to have a low 
complexity ML decoding done over a single partition. We 
note from (8) that any row of the constructed code either 
contains the symbols (s^'s) or its conjugates (with a possible 
sign change). For any h e C Mxl , define a transformation 
denoted by T that takes conjugates of those elements of M x 1 
vector Gm [s]h that contain conjugates of elements of s, and 
we can write 



T{G M [A Mi (s)]h} = £ M ,i(h.)v M ,i( s )> 



(9) 



where £m/$ are M x (M/2) matrices dependent only on 
h, vm/s. are (M/2) x 1 vectors that contain symbols from 
partition i, with i = 1,2. We need a few results from [12] that 
we state here without proof. 



Proposition 1: For any h, s <E C Mxl , 

G m [AmA s )] g m[AmA s )} + 

GI[Am, 2 (s)]Gm[Am,i(s)} =0 m , (10) 

£l A (h)£ Mt2 (h) = M / 2 , 



(11) 



det {G 2 m[A2m,i(s)]} = det {G M [Am, i(sm,i - sm,2)]} 

x dot {G M [-4 M ,i(sAf,i + s M ,2)]}, (12) 

where for any 2M x 1 vector z, we define a trans- 
formation denoted by z that interchanges the two halves 
of z with a sign change for the second half, i.e. z = 
[—zm+i, ■ ■ ■ , —z 2 m, Z\, ■ ■ ■ , Zu\- 

By taking conjugates appropriately, we can derive a mod- 
ified signal model from (3) for receive antenna n, (n = 
1,---,N), as 

X n (s) = [£M.i(H n )v M .i(s) + £M,2{H n )v M ,2{s)] + V n , 

where H n is the nth column of H and X n and V n are derived 
from nth column of X and V respectively by taking the 
conjugates of some or all their elements. Let the singular value 
decomposition (SVD) of £Ms{H n ) be given by 



£M,i{H n ) — UmsSmsW m i , 



(14) 



where Um,i and Wu,i are unitary and Sm,% is a M x (M/2) 
diagonal matrix. Let Sm.i be a M x (M/2) diagonal matrix 
whose diagonal elements are inverse of diagonal elements of 
Sm i and hence 



SmaS 



M,i 



jM/2 qM/2 
qM/2 qM/2 



(15) 



and SmaS m t SM,i — Sma- Multiplying both sides of (13) by 

U M aS M aW m { £ m = U M ,iS M ,iS Mi U Mi , we get 



Uma$>maS m JJ m i X n (s) 



- r £M,i(H n )VM,i(s) 



+U M tl U Ml Vn, (16) 

where we have used £ M l {H n )£M. 2 {H n ) = to cancel the 
contribution of other partition. Note that using (15), it follows 
that V n — UmaSmaS m JLJ m y n has the same the statistics 
as V n - We can rewrite (16) as 



X n (s) 



-,S M aW m t V M ,t(s) + V n 



(17) 



Using (7), one can iteratively generate the equivalent chan- 
nels for each partitions with h M1 = [hi,-- - , h M ] and h M2 = 

[hn+i, 1 ■ ■ j h 2M \, as 



£2M,l(h) = 
^2M,2(h) = 



^M,l(Vl) ^M,2(hM,2) 
^M,l( h M,2) ^M^tVl) 

-£m,2(1im,i) —£m,i(^m, 2 ) 



y M. 2 



(h M ,2) 



^M,l( h M,l) 



(18) 



(19) 



B. Codes with controllable decoding complexity 

Before we get to the code design, we first prove some 
properties that are given in the following propositions. 
Proposition 2: The matrices 

T2M,l(h2M) = ^M.l^M )£2M,l(h2Af), (20) 
T2M,2(h2M) = ^2M,2(' 1 2M)^2M,2(h2M) J (21) 

Kjf(hjlf,l! hM,2) = ^i,/.i(hM,l)^JW,2(hM,2) - 

(h M , 2 )fM, 2 (hM,i), (22) 

YjVf(hM,lj h Mi2 ) = ^M,i(hM,l)^M,l(hM,2) + 

^M ;1 (h M ,2)^M,l(hM,l), (23) 

ZM(hM,l,hM,2) = fl/.2(hM,l)^M,2(hM,2) + 

4/, 2 ( h Af,2)fM,2(h M ,l) (24) 

are real V h 2M G C 2Mxl , h M) i € C Mxl , h M , 2 G C Mxl . 
Proo/: Omitted. ■ 



Proposition 3: For any h.M,i,h-M,2 G C Mxl , 
if Ym(1im,i, h M ,2) and T M] i(hM,i) have the same eigen- 
vectors and Zm0^m,i,^m,2) and Tj,/,2(hM,i) have the same 
eigenvectors, then for any h2M,g2M € C 2Mxl , eigenvectors 
of Y 2 M(h2M, g2M), T 2 M,i(h2M) are the same, and similarly, 
the eigenvectors of Z 2 m0o.2M, S2m), T 2 M,2(h2M) are also 
the same. 

Proof: Omitted. ■ 



Proposition 5: If 



Proposition 4: If for any h^M € C 



4Mxl 



»4M 
b4M 



(25) 



is an eigenvector for T4M,i(h.4M) with A4M as the associated 
eigenvalue, then the eigenvector of T4M,2(h.4M) is 



b4M 

— &4M 

with the same eigenvalue A 4 m- Furthermore, 

^■2M^>4M — A 4M a4M, 



^■2M a 4M — A 4 ^b 4 Af, 



(26) 



(27) 
(28) 



where the dependence of K on the channel realization is 
dropped for convenience. 

Proof: Omitted. ■ 



&4M 
b4M 



(29) 



is an eigenvector of T4M,i(h4Af), then the eigenvectors of 
T 8 M,i( h 8M) are 



1 

71 



&4M 




&4M 


b4M 


and — ^ 


b4M 


b4M 




— b 4 M 


— &4M 




&4M 



(30) 



Proof: We prove this by induction. It is easy to check it 
for M = 8. Let us assume that this is true for T^ ; i(hfe) V k < 

4M i.e. if ^ 4M is an eigenvector of Tam i(h4M), then 
a 4M = f 2M is an eigenvector of T 2 m,i (h 2 M). By using 

t>2M 

Proposition 2, T 8 M,i(hsM) can be written as into smaller 
parts as 

TsM,l(h8M) = 



T4A// j i(h4M,l) + T4M,l(h4M,2) 

K4M(h4M,l, h4M,2) 
K 4 r M (h4M,l, h4M,2) 

T4Af i 2(h4M,2) + T4M i 2(ll4M,l) 



(31) 



We have to show that if 
T 4M ,i(h 4 M), then 



T8M,l(h-8M) 



From the induction assumption, 



&4M 
b4M 



&4M 




&4M 


b4M 


— Asm 


b4M 


b4M 


b4M 


— a4M 




— a 4 M 



&4M 
b4Af 



is an eigenvector of 



(32) 



is an eigenvector 



of T 4M ,i(h4M,i) and T 4 M,i(h4M,2) with eigenvalues X iM , 
Af M respectively and using Proposition 4, these are also the 
eigenvectors of T 4M ,2(h4M,i) and T 4 M,2(h4M,2)- Substitut- 
ing in (31) and (32), we have to show that 



K.4M (h-4M, 1 , h4M,2 ) 



If 



K 4M (h.4M,l, h4M,2) 

&4M 
b4A^f 



b4M 

— a4M 

&4M 
b4M 



A 8M 



— A 8M 



»4M 
b4M 

b4M 
-a4M 



(33) 
(34) 



is an eigenvector of T4M,i(h4M), then it follows 



from the induction assumption that s^m is an eigenvector of 
T2M,i (and Y2M) and b4M is an eigenvector of T2M,2 (and 
Z2Af), where the dependence of T on the channel realization 
is dropped for convenience. Using Propositions 3 and 4, we 
have 

K2M (h2M,l, h-2M,2)b4M = A 4M (b-2M,l , h-2M,2)a4M , (35) 
K 2M (h2M,l, h 2 M,2)a4M = A4 M (h 2 M4, h 2 M,2)b4M , (36) 
Y2Af (h 2 M,l, h 2 Af,2)a4M = A2 M (h 2 M,l, h 2 M,2)a4M, (37) 
Z2M(h-2M,li h 2 M,2)a4M = A 4M (h 2 M,l , h-2M,2)b4M (38) 



Note that ^ K 4M (h 4¥A , h 4M>2 ) 

-K2M(h2M,l, tl2M,3) + K2M(h2M,4- h-2M,2) 

— Y2M (h2M,l ) h 2 M,4) + (h 2 M,2; h2Af ,3) 

— Z2Af (h2M,2 , h2M,3) + Z 2M (h2M,l, h2M,4) 

Kj M (h2M,l, h2M,3) — (h 2 M,4 5 h 2 M,2) 

using (35), (36), (37), (38), we have 



, hence 



K4M(h4M,l: h4M,2) 



b4M 



— A x 



8M 



b4M 



, (39) 



where Ag M — — A4 M (h 2 M,i, h 2 M,3) + A4 M (h 2 M,4, h 2 M,2) + 
A4 M (h 2M ,i,h 2M ,4) - A| M (h 2M ,2,h 2M ,3)- This proves (33). 
Hence 



-K2M(h2M,l, h2M,3) + K2A//(h2M,4, h2A'/,2) 

-Y2M (h2M,l, h2M,4) + (h2M,2, h2M,3) 

— Z2M(h2M,2j h2M,3) + Z2 M (h2M,l, h2M,< 



K 2M (h: 



2M,1 , n2A/,3 J 



K2 M (h 2 M,4, h 2 



— A 8M 



b4M 
— S4M 

b4M 




By interchanging the first half of the rows with the second 
half, then interchanging the first half of the columns with the 
second half, then multiplying the first half of columns and the 
second half of rows with —1, and using the fact that Ym, Zm 
are real, Hermitian matrices, we can write the above equations 
as 



-Kj M (h 2 M,l , h 2 jy/,3) + K* M (h 2 M,4, h 2 M,2) 



iJVJ . ^~ , ■ 

-Y2 M (h2M,l, h2M,4) + Y2 M (h2M,2: h2M,3) 

K2M (h2M,l, h2M t 3) — K2M(h2M,4, h2M,2) 



— A 8M 



&4M 
b4M 

b4M 

-&4M 



or 



K4 M (h4M,l, h4M,2) 



MM 

Hence if 

&4M 
b4M 
b4M 

-&4M 



»4M 
b4M 



A 8M 



b4M 
-S4M 



(40) 



&4M 
b4M 



is an eigenvector of T4A-/;i(h4M), 
is an eigenvector of T 8 m,i0^8m)- Similarly, by 



using (33) and (34), one can show that 
eigenvector for T 8 m,i(^8m)- Q-E.D. 



»4Af 
b4M 

»4Af 



is also an 



Example: For M = 4, the eigenvector matrix for T4.1 (or 
lputec 

W 4A = 



Wm,i in (14)) is computed as 

1 

and using Proposition 5, the eigenvector matrix for T$.i is 
given by 



(41) 



W 8 ,i = 



1 



1 1 1 

1 -1 -1 

-1 -1 1 

1 -1 1 



(42) 



Decoding complexity 



Fig. 1. Diversity versus decoding complexity tradeoff for the proposed rate 
1 codes. 



An important aspect of the eigenvector matrices of Tm,i 
and Tm,2 is that they are independent of the channel realiza- 
tion. This property is quite useful in constructing codes with 
controllable ML decoding complexity. 

Proposition 6: The eigenvalues for the first partition of the 
code matrix are given by 

vmA s *)} (43) 

where Qm.\ = \/M/2Wm,i and the determinant for the first 
partition of the code matrix is given by 

det {G m [Am,i(s)}} = f (<2L,i«m,i(s)) (44) 

and for any n length vector q, /(q) = Ylk=i \li\ 2 - 

Proof: Omitted. ■ 



Result is similar for the second partition, and due to similarity 
with the above Proposition, we omit it. 

We note here from (17) that it is Wm,% that dictates the 
ML decoding complexity. For example, we could precode the 
information-carrying symbol vector c in (1) such that 



%,i(c) = W M ,iVM,i(s) 



(45) 



According to (17), this code will admit single symbol ML 
decoding. But this will give the determinant of the error code 
matrix as 

M 

det {GmIAmAWm Ac - e))}} = — f(v M ^c-e)) (46) 

Since the elements of c and e are drawn from the same 
constellation, hence even if any element of c and e is the 
same, then 

min det {Gm[Am.i(Wm Ac -e))}} = (47) 

c.e.c^e 

For single symbol decoding, the minimum rank of Gm would 
be 2 (also the rank of G^ m Gm)- In general, ML decoding of 



Proposed code 
^ — Ideal orthogonal code 
No fading 




O 2 symbols 
O 4 symbols 
8 symbols 

Ideal orthogonal code 




12 14 16 
SNR (dB) 



22 



Fig. 2. SER versus SNR for various M and N = 1 with QPSK modulation 
for the rate-1 constructed codes with M/2 symbols decoded together and the 
ideal orthogonal codes. 



Fig. 3. SER versus SNR for M = 16 and N = 1 with QPSK modulation for 
the rate-1 constructed and the ideal codes for varying decoding complexity. 



M/2™, n — 1, 2, • • • , log 2 M, symbols together would mean 
having the precoding matrix as a block diagonal matrix with 
each block as Wm/2"- 1 ,! (scaled appropriately) with constel- 
lation rotation to ensure that the rank of Gm[Am,i( c — e )] is 
M/2"- 1 . 

The only way to achieve full diversity would be to choose 
n = 1 or decode M/2 information symbols together. One can 
employ various methods like constellation rotation given in 
[9], [11], [14], [15]. 

By using this block diagonal structure, we can construct 
codes with ML decoding of different symbols together and 
hence the ML decoding complexity can be controlled. We plot 
the diversity versus complexity tradeoff in Fig. 1, where M 
is not necessarily a power of 2. The code design for such M 
is done by consturcting a code for 2 ^ l ° S2 M l transmit antennas 
that admits ML decoding complexity of 2 r io S2 A/ l -« ( n — 
1, 2, • • •) and has rank of each partition as 2 r io S2 Ml-n+i an( j 
then deleting columns suitably chosen to retain the same rank 
and to have the code matrix of the size 2 l" log2 M l x M. 

Note that it is not necessary to assume that P = M i.e. 
unit rate. One could design codes with P < M that may have 
additional coding gain while sacrificing code rate. 

IV. Numerical Results 

The symbol error rate (SER) versus the average SNR per 
receive antenna for the proposed rate-1 code that admits 
decoding in pairs of symbols is plotted in Fig. 2 with QPSK 
modulation for M = 4,8, 16 and N = 1. Also plotted is the 
performance of an ideal rate-1 orthogonal space-time codes 
(non-existent for M > 2) with equivalent channel as . 

Fig. 3 plots the SER curves for M — 16 and N = 1 for 
different ML decoding complexities. Fig. 1 plots the diversity 
versus complexity tradeoff for the proposed codes. 

V. Conclusions 

We have constructed a class of linear space-time codes that 
have controllable ML decoding complexity for any number of 
transmit antennas. The diversity versus decoding complexity 
tradeoff is shown. We show that one can design rate 1 



codes that achieve performance quite close to the rate 1 ideal 
orthogonal codes (non-existent for for M > 2). 
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